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Abstract 

In this work we study the phenomenological consequences of the environment dependence of 
neutrino mass on solar and reactor neutrino phenomenology. We concentrate on mass varying 
neutrino scenarios in which the enviroment dependence is induced by Yukawa interactions of a 
light neutral scalar particle which couples to neutrinos and matter. Under the assumption of 
one mass scale dominance, we perform a global analysis of solar and KamLAND neutrino data 
which depends on 4 parameters: the two standard oscillation parameters, Am 2 , 2 ± and tan 2 9±2, and 
two new coefficients which parameterize the environment dependence of the neutrino mass. We 
find that, generically, the inclusion of the environment dependent terms does not lead to a very 
statistically significant improvement on the description of the data in the most favoured MSW LMA 
(or LMA-I) region. It does, however, substantially improve the fit in the high-Am 2 LMA (or LMA- 
II) region which can be allowed at 98.9% CL. Conversely, the analysis allow us to place stringent 
constraints on the size of the environment dependence terms which can be translated on a bound 



on the product of the effective neutrino-scalar (A") and matter-scalar (A ) Yukawa couplings, as 
a function of the scalar field mass (m s ) in these models, \\ v \ N \ ( 10 ~ 7 s eV ) < 3.0 x 1CT 28 (at 90% 
CL) . 
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I. INTRODUCTION 



The possibility of environment dependence (ED) of the effective neutrino mass was first 
proposed as a possible solution to the solar neutrino deficit by Wolfenstein . In the Stan- 
dard Model, the vector part of the standard charged current and neutral current neutrino- 
matter interactions contribute to the neutrino evolution equation as an energy independent 
potential term. The potential is proportional to the electron density and it has different 
sign for neutrinos and anti-neutrinos, even in the CP conserving case. It gives rise to the 
well-known Mikheyev-Smirnov- Wolfenstein (MSW) effect 0, Q| which is crucial in the inter- 
pretation of the solar neutrino data. 

In most neutrino mass models, new sources of ED of the effective neutrino mass arise as a 
natural feature due to the presence of nonstandard neutrino interactions with matter & If 
the new interaction can be cast as a neutral or charged vector current, it will also contribute 
as an energy independent potential to the neutrino evolution equation. The phenomenolog- 
ical implications of such non-standard interactions in neutrino oscillations have been widely 
considered in the literature nun 

New physics in the form of Yukawa interactions of neutrinos and matter with a neutral 
light scalar particle modify the kinetic part of the neutrino evolution equation. Such neutral 
scalar interactions induce a dependence of the neutrino mass on the environment 0| which 
is energy independent and has the same sign for neutrinos and anti-neutrinos. Recently, this 
form of ED of the neutrino mass has recefved renewed attention after Ref. | discussed the 
possibility that such mass varying neutrinos (MaVaNs) can behave as a negative pressure 
fluid which contributes to the origin of the cosmic acceleration. This scenario establishes 
a connection between neutrino mass and dark energy with interesting cosmological conse- 
quences 0, 0, E3 

In the MaVaNs scheme presented in Ref. , the neutrino mass arises from the interaction 
with a scalar field, the acceleron, whose effective potential chang function of the 

neutrino density. As a consequence, the neutrino mass depends on the neutrino density in 
the medium. A subsequent work, Ref. , investigated the possibility that neutrino masses 
depend on the visible matter density as well. Such a dependence would be induced by non- 
renormalizable operators which would couple the acceleron to the visible matter. This form 

so lead to interesting phenomenological consequences for 



of ED of the neutrino mass could a 
neutrino oscillations 
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For solar neutrinos, in Ref. 3] it was shown that, generically, due to the dependence of 
the neutrino mass on the neutrino density, these scenarios establish a connection between 
the effective Am 2 in the Sun and the absolute neutrino mass scale. Due to this effect, 
the description of solar neutrino data worsens for neutrinos with degenerate masses. On 
the other hand, for hierarchical neutrino masses the dominant effect is the dependence of 
the neutrino mass on the visible matter density. In Ref. 0] it was shown that for some 
particular values of the scalar-matter couplings this effect can improve the agreement with 
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solar neutrino data. 

In this article we investigate the characteristic effects of the dependence of the neutrino 
mass on the matter density for solar neutrinos and reactor antineutrinos. We perform a 



15. 
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and KamLAND data 



2l| in these scenarios. 



combined analysis of solar 
Our results show that: (i) the inclusion of the ED terms can lead to certain improvement of 
the quality of the fit in the most favour LMA-I re gion for well determined values of the new 
parameters, (in agreement with the result of Ref. |l3|), but this improvement does not hold 
much statistical significance; (ii) the inclusion of these effects, can substantially improve the 
quality of the fit in the high-Am 2 (LMA-II) region which can be allowed at 98.9% CL; (Hi) 
generically, the combined analysis of solar and KamLAND data results into a constraint on 
the possible dependence of the neutrino mass on the ordinary matter density. 

In Sec. |nj we introduce the general theoretical framework which we will consider in this 
paper. In Sec. IIIIl we discuss how the ED modifies neutrino oscillations in matter. In 
Sec. I1V1 we examine how these modifications can affect the current allowed solar neutrino 
oscillation parameter region and establish the constraints that solar and reactor neutrino 
data can impose on the couplings to the scalar field. Finally, in Sec. El we discuss our results 
and summarize our conclusions. 



II. FORMALISM 

For the sake of concreteness, we consider here an effective low energy model containing 
the Standard Model particles plus a light scalar (<fi) of mass mj which couples very weakly 
both to neutrinos (i>j) and the matter fields f — e, n,p. 

The Lagrangian takes the form 

£ = J> - m°) + f W - m °f) f+\ [<!> (° 2 ~ m2 s) 4\ +E Xf ff^ » 

* / i'j f 

(1) 

where mj are the vacuum mass that the neutrinos would have in the presence of the cosmic 
neutrino background. A^- and A* are, respectively, the effective neutrino-scalar and matter- 
scalar couplings. We have written a Lagrangian for Dirac neutrinos but equivalently it could 
be written for Majorana neutrinos. 

In a medium with some additional neutrino background (either relativistic or non- 
relativistic) as well as non-relativistic matter (electrons, protons and neutrons), neutrinos 
acquire masses which obey the following set of integral equations 

mij(r) = m°Sij - Mij(r), 
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nf(r) is the number density for the fermion /, and f a (r,k) is the sum of neutrino and 
antineutrino "a" occupation numbers for momentum k in addition to the cosmic background 
neutrinos. 

In the context of the dark energy-related MaVaNs models of Ref the scalar (j) would 

be the acceleron - with mass in the range ms ~ 10~ 6 -10~ 8 eV - which, when acquiring a 
non- vanishing expectation value, (0), gives a contribution to the neutrino mass. This in 
turn implies that the acceleron effective potential receives a contribution which changes as 
a function of the neutrino density, so that 

^ _ drn^ 



(3) 

(*) 

\', the effective low energy couplings of the acceleron to visible matter, come from non- 
renormalizable operators which couple the acceleron to the visible matter, such as might 
arise from quantum gravity. They can be parametrized as 



where Mm is the Plank scale. Tests of the gravitational inverse square law require X n , X p 



< 



10 



-21 



22j for any scalar with ms > 10 -11 eV. 
The results in Eq.(J2J) and Eq.Q correspond to the first order term in the Taylor expansion 
around the present epoch background value of <p. In general, for the required flat potentials 
in these models, one needs to go beyond first order and the neutrino mass is not linearly 
proportional to the number density of the particles in the background. The exact dependence 
on the background densities is function of the specific form assumed for the scalar potential. 
This is mostly relevant for the neutrino density contribution to the neutrino mass. 

It has been recently argued j^] that, generically, these models contain a catastrophic 
instability which occurs when neutrinos become non-relativistic. As a consequence the 
acceleron coupled neutrinos must be extremely light. Thus, in what follows we assume the 
vacuum neutrino masses to be hierarchical 

= m ° < m ° < m° . (5) 

For solar neutrinos of hierarchical masses, as discussed in Ref. ["ill Q], the dominant 
contribution to the neutrino mass is due to the matter background density. Correspondingly, 
we neglect the contribution to the neutrino mass from the background neutrino density and 
we concentrate on the matter density dependence: 

M *^) = JrE A SM- (6) 



s 



f 

This is very similar to the scenario considered in Ref. \^ . 
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Finally, let's mention that by assuming Eq.(^) we do not consider the possibility of 
additional light mixed sterile neutrinos which may appear in some specific realizations of 
MaVaNs scenarios 24J and which can lead to other interesting effects in oscillation neutrino 
phenomenology and cosmology 24 , 25 , 2f 



III. EFFECTS IN SOLAR NEUTRINO OSCILLATIONS 



Q, K2K Q, solar Q, 
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The minimum joint description of atmospheric [2 
20l | and reactor |2jJ, |29j data requires that all the three known neutrinos take part in the 
oscillations. The mixing parameters are encoded in the 3x3 lepton mixing matrix which 
can be conveniently parametrized in the standard form 



where c,- 
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According to the current data, the neutrino mass squared differences can be chosen so 
that 



Am| 



Am^ < \Aml 



31 1 



|Am^ 2 | 



Am 



atm ' 



(8) 



As a consequence of the fact that Am^/IAm^l w 0.03, for solar and KamLAND neutrinos, 
the oscillations with the atmospheric oscillation length are completely averaged and the 
interpretation of these data in the neutrino oscillation framework depends mostly on Am^, 
6 12 and 6*13, while atmospheric and K2K neutrinos oscillations are controlled by A mji, 
6*23 and #13. Furthermore, the negative results from the CHOOZ reactor experiment [29j 
imply that the mixing angle connecting the solar and atmospheric oscillation channels, #13, 
is severely constrained (sin 2 $13 < 0.041 at 3a Altogether, it is found that the 3-u 

oscillations effectively factorize into 2-v oscillations of the two different subsystems: solar 
and atmospheric. 

With the inclusion of the ED terms (Eq. (JHJ)) it is not warranted that such factorization 
holds. We will assume that this is still the case and study their effect on solar and KamLAND 
oscillations under the hypothesis of one mass-scale dominance. Under this assumption, we 
parametrize the evolution equation as fl^ : 



dr 



2E U 



U ei2 



Ml(r) 
M|(r) 



Mf(r) 




(9) 



where we have assumed the neutrinos to follow the hierarchy given in Eq.flHJ). Here Vcc( r ) — 
V2G F n e { r) is the MSW potential proportional to the electron number density n e (r) in the 
medium. Ug 12 is the 2x2 mixing matrix in vacuum parameterized by the angle 612, and 
Mj(r) are the ED contributions to the neutrino masses. 
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In general, fo r gi ven matter density profiles, Eq. Q has to be solved numerically. As 
discussed in Ref. [l3| in most of the parameter space allowed by KamLAND and solar data, 
for all practical purposes, the transition is adiabatic and the evolution equation can be solved 
analytically to give the survival probability 

P ee = ^ + icos20™ 12 cos20 12 , (10) 

where #™ 12 is the effective neutrino mixing angle at the neutrino production point r in the 
medium, explicitly given by 

ofj m (AM|(tq) cos2go^ - 2E u V cc (r )) 



(AMl(r ) cos20 o ,i 2 - 2E u V cc (r )y + (AMf^ro) sin20 o ,i2) 2 



with 



AM 2 >o) = M 3 4 (r ) + (^MM^ (12) 



cos 202i - M£(r ) sin 20i2 



cos 20 o .i2 = i (13) 

W^fM) 2 



and where 



AJW*(r ) = (m° 2 - M 2 (r )) 2 - M 2 (r ). (14) 

As discussed in the previous section, in general, Mj(r) can be an arbitrary function of the 
background matter density. For the sake of concreteness we will assume a dependence in ac- 
cordance with the results obtained in the linear approximation given in Eq.©. Furthermore, 
using Eq.(0J), A e <C \ n = \ p = \ N , so we will parametrize these terms as: 

p{r) 



Mi(r) = <*i 



(gr/cm 3 ) 



(15) 



where p is the matter density, and from Eq.© we find the characteristic value of the a 
coefficients to be 2 

a~4.8xl0 23 A v A JV (^1) e V . (16) 

V m s J 

One must notice, however, that, as long as the transition is adiabatic, the survival prob- 
ability only depends on the value of Mj(r) at the neutrino production point. Therefore it 
only depends on the exact functional form of Mj(r) via the averaging over the neutrino 
production point distributions. 

The survival probability for anti-neutrinos, Peg, which is relevant for KamLAND, takes 
the form 

P ge - = l-sin 2 20 o ": i2 sin 2 (^|^) , (17) 
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where cos 2#™ 12 is defined as in Eq. (II lj) and Am 2 KL is the denominator of this equation but 
replacing Vcc by —Vcc and assuming a constant matter density p ~ 3 gr/cm 3 , typical of 
the Earth's crust. 

To illustrate the effects of the a coefficients, we show in Fig. ^ the evolution of the mass 
eigenvalues ni\ and m<i in matter as a function of VqqE v for different values of (keeping 
ol\ = «3 = 0). As a reference, we also show in this figure the standard MSW evolution 
curve (solid line) for the oscillation parameters at A?72q 21 — 

( m 0) 2 = 8 x 10~ 5 eV 2 and 

tan 2 9 12 = 0.4, a point which explains very well both solar and KamLAND data. From this 
plot we can appreciate that in the region relevant to solar neutrino experiments the evolution 
of the mass eigenvalues is not significantly different from the MSW one if |gk2 j ^ 10~ 5 eV. 
For larger values of «2, such as jo^l = 10~ 4 eV, we expect solar neutrinos to be affected. On 
the other hand, KamLAND data is very little affected by the ED terms in this range of ot^. 

Figure 1 also illustrates a curious feature of these scenarios: the fact that it is possible to 
find a value of the matter dependence term which exactly cancels Am§ 21 . It can be seen, 
directly from Eqs. 0, that if for a particular point, r , in the medium, m 2 = M 2 (ro) and 
Mi = M 3 = (ai = 0:3 = 0) the lower mass eigenstate will be zero while the higher one will 
be at the corresponding value of 2Vcc( r o)E u . 

Non-adiabatic effects in the Sun can also occur. In the region of relatively small a 
parameters, non-adiabaticity occurs when the parameters are "tunned" to give a vanishing 
effective Am 21 (the denominator of Eq. (JUJ). This can be achieved, for example, with 
a± = a.2 = by solving the following set of equations inside the Sun: 



It can be shown that for = z5.5 x 10 -5 eV, tan 2 9 = 0.3 and E u = 10 MeV this set 
of equations are fulfilled at r/R & ~ 0.027, and the neutrinos would suffer a non-adiabatic 
transition on their way out of the Sun. However, in general for the small values of the a 
parameters discussed here, these non-adiabatic effects do not lead to a better description of 
the solar neutrino data. 

More generically, non-adiabatic effects occur for sufficiently large values of the a param- 
eters so that one can disregard the standard MSW potential Vcc an d the vacuum mass 
m 2 with respect to the matter density mass dependent terms. In this case, as seen from 
Eq. (jllj) . the mixing angle inside the Sun is constant and controlled by the a's. At the 
border of the Sun, as the density goes to zero, the mixing angle is driven to its vacuum value 
in a strongly non-adiabatic transition. This scenario would be equivalent to a vacuum-like 
oscillation for solar neutrinos with the ED of neutrino mass having to play a leading role in 
the interpretation of terrestrial neutrino experiments. We will leave the detailed analysis of 
the consequences of this type of non-adiabatic transitions for a future work. 



(m 2 , ) 2 cos2#-2M 3 2 (r)sin2# = 2E u V cc (r), 
(m 2 ) ) 2 sin2# + 2M 3 2 (r)cos2# = 0. 



(18) 
(19) 
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FIG. 1: Evolution of the neutrino mass eigenvalues in matter. The solid lines represent the familiar 
MSW evolution. In this plot we have fixed (to!]) 2 = 8 x 1CT 5 eV 2 , tan 2 6*12 = 0.4 and ot\ = 03 = 0. 
In the upper (lower) panel the dotted line represents the case Q2 = — 10 -5 eV (02 = — 10 -4 eV) 
and the dashed line the case «2 = +10 -5 eV (02 = +10 -4 ) eV. The shaded regions correspond 
to typical values of Vcc in neutrino production region in the center of the Sun for the Solar z/'s 
region, and a constant Earth crust density of 3 g/cm 3 , with a proton density fraction oiY = 0.5 
and neutrino energies varying from 3 to 10 MeV for the KamLAND region. 

IV. CONSTRAINTS FROM SOLAR AND REACTOR NEUTRINO DATA 

We present in this section the results of the global analysis of solar and KamLAND for 
the specific realization discussed in the previous section. Furthermore, for simplicity, we will 
restrict ourselves to the case Mi(r) = mi = 0. 

Details of our solar neutrino analysis have been described in previous papers 
use the solar fluxes from Bahcall and Serenelli (2005) j^. The solar neutrino data includes 



total of 119 data points: the Gallium [lfi[ll7 



321. We 

a 

and Chlorine [l5j (1 data point) radiochemical 
rates, the Super-Kamiokande |18( zenith spectrum (44 bins), and SNO data reported for 
phase 1 and phase 2. The SNO data used consists of the total day-night spectrum measured 
in the pure D 2 (SNO-I) phase (34 data points) [3], plus the full data set corresponding 
to the Salt Phase (SNO-II) [2^. This last one includes the NC and ES event rates during 
the day and during the night (4 data points), and the CC day-night spectral data (34 
data points). The analysis of the full data set of SNO-II is new to this work. It is done 
by a x 2 analysis using the experimental systematic and statistical uncertainties and their 
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correlations presented in |20j , together with the theoretical uncertainties. In combining with 
the SNO-I data, only the theoretical uncertainties are assumed to be correlated between the 
two phases. The experimental systematics errors are considered to be uncorrelated between 
both phases. 

For KamLAND, we directly adapt the y 2 map as given by the KamLAND collaboration 
for their unbinned rate+shape analysis [34j which uses 258 observed neutrino candidate 
events and gives, for the standard oscillation analysis, a The corresponding 

Baker-Cousins x 2 f° r the 13 energy bin analysis is Xmm = 13.1/11 dof. The effect of MaVaN's 
parameters in KamLAND result was calculated assuming a constant Earth density of 3 
g/cm 3 , and assuming that KamLAND are sensitive to the vacuum value of Awiq 21 and #i 2 
through an effective mass and mixing in a constant Earth density, respectivelly given by the 
denominator of Eq. (|TT|) and Eq. (fTB^) . as described in Eq. (fTTj) . 

In presence of the ED contribution to the masses, the analysis of solar and KamLAND 
data depends on four parameters: the two standard oscillation parameters AwIq 21 = (m^) 2 , 
and tan 2 # 12 , and the two ED coefficients, a 2 , and a 3 . In this case, in order to cover the full 
CP conserving parameter space we allow the a parameters to vary in the range 

—00 < a 2 < 00 — 00 < a 2 < 00 (20) 

We find the best fit point 

tan 2 On = 0.49 Am 2 21 = 8.4 x 10" 5 eV 2 
a 2 = 10" 4 eV a 3 = % 2.0 x 10~ 5 eV . (21) 

This is to be compared with the best fit point for no ED of the neutrino mass, i.e., 
a 2 = ce 3 = 

tan 2 12 = 0.44 Am 2 21 = 7.9 x 10" 5 eV 2 . 

A X 2 = 2.5 (22) 

where A% 2 is given with respect to the minimum at the best fit point in Eq. (J21J) . Thus 
we find that although the inclusion of the ED terms can lead to a small improvement of 
the quality of the fit (in agreement with the result of Ref. 0]), this improvement is not 
statistically very significant leading only to a decrease of A% 2 = 2.5 even at the cost of 
introducing two new parameters. 

We show in Fig. |2] the result of the global analysis of solar data plus KamLAND data 
in the form of the allowed two-dimensional regions at 3cr CL in the (Am 21 , tan 2 # 12 ) plane 
after marginalization over a 2 and a 3 . The standard MSW allowed region is also showed for 
reference. As seen in the figure, allowing for ED of the neutrino masses enlarges only slightly 
the allowed range of Am^ and tan 2 # 12 in the LMA-I region. In contrast to the standard 
MSW analysis, where the limits on the mixing angle come basically from solar neutrinos, 
here it is KamLAND data that control the lower limits for the mixing angle. 
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FIG. 2: Allowed regions from the global analysis of solar plus KamLAND data in the 
(Amp tan 2 612) parameter space at 3a CL (2dof). The best fit point at tan 2 9\2 = 0.5 and 
Attt-q 21 = 8.7 x 10~ 5 eV 2 is represented by the star. The standard MSW allowed region is also 
shown for reference. 



Most interestingly, we also find that the description of the solar data in the high-Am 2 
(LMA-II) region can be significantly improved so there is a new allowed solution at the 
98.9% CL. The best fit point in this region is obtained for 

tan 2 0i 2 = 0.5 Am 2 21 = 1.75 x 10" 4 eV 2 
a 2 = 1.3 x 10" 4 eV a 3 = i 2.0 x 10" 5 eV 

A X 2 = 8.9 (23) 

While this region is excluded at more than 4cr for standard MSW oscillations, it is allowed 
at 98.9% CL (2.55a) in the presence of environmental effects with | CK3 1 < 3.2 x 10~ 5 and 
2.8 x 10~ 5 < «2 < 2.0 x 10~ 4 . Basically the CL at which this region is presently allowed 
is determined by KamLAND data 2l( because the fit to the solar data cannot discriminate 
between the LMA-I and LMA-II regions once the ED terms are included. Clearly this implies 
that this solution will be further tested by a more precise determination of the antineutrino 
spectrum in KamLAND. 

We show in FigOlthe survival probability for this best fit point in the high-Am 2 (LMA-II) 
region in the presence of ED effects together with the extracted average survival probabilities 
for the low energy (pp) , intermediate energy ( 7 Be, pep and CNO) and high energy solar 
neutrinos ( 8 B and hep) from Ref. 0|. For comparison we also show the survival probability 
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FIG. 3: f e survival probability in the Sun versus neutrino energy for the best fit point in the high- 
Am 2 region in the presence of ED effects (Eq. (|23|)). The dotted line is the survival probability 
for conventional oscillations (c^ = 0) with the same values of Am^g and #12. These survival 
probabilities have been obtained for neutrinos produced around x = 0.05 as it is characteristic of 
8 B and 7 Be neutrinos. The data points are the extracted average survival probabilities for the low 
energy (pp) , intermediate energy ( 7 Be, pep and CNO) and high energy solar neutrinos ( 8 B and 
hep) from Ref. [x3 | . 

for conventional oscillations (c^ = 0) with the same values of Amfio ano ^ ^12- From the 
figure it is clear that the inclusion of the ED parameters, leads to an improvement on the 
description of the solar data for all the energies being this more significant for intermediate- 
and high-energy neutrinos. 

On the contrary, unlike for the case of non-standard neutrino interactions discussed in 
Ref. 0,0, the low- Am 2 (LMA-0) region is still disfavoured at more than 3 sigma by the 
global KamLAND and solar data analysis even in the presence of the new "kinetic-like" ED 
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FIG. 4: Allowed regions from the global analysis of solar and solar plus KamLAND data in the 
(02,03) parameter space. The curves correspond to la, 90%, 95%, 99% and 3a CL (2dof). The 
best fit point at 02 = 10~ 4 eV and 03 = i 2.0 x 10~ 5 eV, represented by a star, is also shown. 



effects discussed here. This is due to the different energy dependence of the new physics 
effects in the two cases. For the "kinetic-like" effects it is not possible to suppress matter 
effects in the Earth for the high energy neutrinos (to fit the SK and SNO negative results 
on the day-night asymmetry within LMA-0) without spoiling the agreement of the survival 
probability at intermediate energies with the result of the radiochemical experiments. 

Conversely, the global analysis of solar and KamLAND data results into the constraint 
of the possible size of the ED contribution to the neutrino mass. This is illustrated in Fig0] 
where we show the result of the global analysis in the form of the allowed two-dimensional 
regions in the (a2, 03) parameter space after marginalization over Amg 21 , tan 2 #12. The full 
regions correspond to la, 95% and 3 a CL while the curves correspond to 90 and 99% CL. 
As seen in the figure, for CL > l.ler the regions are connected to the 02 = 03 = case 
and they are always bounded. In other words, the analysis show no evidence of any ED 
contribution to the neutrino mass and there is an upper bound on the absolute values of the 
corresponding coefficients. 

Our previous discussion at the end of Sec. IIHI on the behavior of the mass eigenvalues 
shown in Fig.^when M2 — > m®, is the reason behind the exclusion of the narrow gulf around 
«2 ~ 10~ 4 eV. In this region, unless 03 is negative and not too small, it is not possible to 
explain solar neutrino data. 
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In order to quantify the bound on MaVaN's parameters, we display in FigJS] the depen- 
dence of the global x 2 on a 2 (0^3) after marginalization over Am§ 21 , tan 2 # 12 an d a 3 (^2)- 
From the figure we read the following 90% CL (3 a), bounds (with ldof) 

- 2.2 x 1(T 5 < a 2 /eV < 1.4 x 10~ 4 (-5.6 x 10" 5 < a 2 /eV < 1.7 x 10" 4 ) (24) 
H/eV < 2.3 x 10~ 5 (Kl/eV < 8.4 x 10" 5 ) for > (25) 
|a 3 |/eV < 3.4 x 10~ 5 (|ar 3 ]/eV < 5.2 x 10" 5 ) for aj < (26) 

These bounds can be converted into a limit on the product of the characteristic effective 
neutrino-scalar and matter-scalar couplings. For example, at 90% CL, 

I a- A"| ( ^—^.) 2 < 3.0 x 10- 28 . (27) 



m s 

We can compare this bound with those derived from tests of the gravitational inverse square 
law (ISL) which require the coupling of the scalar to nucleons |A | < 10~ 21 22] for any 
scalar with ms > 10~ n eV. Thus we find that if the scalar also couples to neutrinos with 
coupling 

A'>3.0xl0- 7 (-^-) 2 (28) 



.10" 7 eV, 

the analysis of solar and KamLAND data yields a more restrictive constraint on the matter- 
scalar couplings than ISL tests. 

Finally, we want to comment on the possible mo del- dependence of these results. There 
are two main sources of arbitrariness in our derivations: the choice of mi = Mi(r) = 
and the assumption that the are linearly dependent on the matter density. Indeed their 
effect is the same: departing of any of these assumptions results into a different functional 
dependence of the effective neutrino masses with the point along the neutrino trajectory, 
rriiir). 

As discussed in the previous section, the basic assumption behind our results is that 
neutrino evolution in matter is adiabatic. As long as this is the case, the survival probabil- 
ity only depends on the value of the effective neutrino masses at the neutrino production 
point and the final results depend very mildly on the exact functional form of m^r). As a 
consequence the generic results will still be valid: there will be a slight improvement on the 
quality of the fit in the LMA-I region, there will be a substantial improvement of the quality 
of the fit on the LMA-II region and generically the combined analysis of solar and Kam- 
LAND data will result into a bound on the strength of the new contributions. Of course, the 
exact numerical values of the corresponding ED couplings will be different. But the order 
of magnitude of the bound on the product of the Yukawa couplings in Eq. pTj) will hold. 



V. DISCUSSION 

We have investigated the phenomenological consequences of a scalar induced ED of the 
effective neutrino mass in the interpretation of solar and reactor neutrino data. For the 
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lex I la I 

FIG. 5: Dependence of x 2 on the ED parameters «2 (upper panel) and 03 (lower panel) after 
marginalization over the other three parameters. 
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sake of concreteness, we consider an effective low energy model containing the Standard 
Model particles plus a light neutral scalar (<ft) of mass mj which couples very weakly both 
to neutrinos (z/j) and the matter fields / = e,n,p. This is described in Sec. HP and its 
consequences to neutrino oscillations in the Sun is discussed in Sec. IIHI 

Assuming the neutrino masses to follow the hierarchy = m\ < m\ < mjj, we have 
performed a combined analysis of the solar neutrino data (118 data points) and KamLAND 
(17 data points) in the context of this effective model. Our analysis, which is described in 
Sec. IIV1 depends on 4 parameters: the two standard oscillation parameters AmJ 21 = (m 2 ) 2 , 
and tan 2 # 12 , and the two ED coefficients, an d «3- We found the best fit point at: 
tan 2 #i2 = 0.49, Am 2 21 = 8.4 x lCT 5 eV 2 , a 2 = 10" 4 eV and a 3 = j2.0x 1(T 5 eV. This 
point corresponds to a decrease of Ax^^ = —2.5 in comparison to the minimum in the 
case where no ED is considered. We conclude that in spite of the inclusion the two extra 
parameters, the improvement of the quality of the fit in the most favoured LMA-I MSW 
region is not very statistically significant. 

Most interestingly, we find that the description of the solar data in the high-Am 2 (LMA- 
II) region can be significantly improved and there is a new allowed solution at the 98.9% CL. 
The best fit point in this region is obtained for tan 2 #12 = 0.5, Am 2 , 21 = 1.75 x 10~ 4 eV 2 , 
a 2 = 1-3 x 10 -4 eV and = i 2.0 x 10 -5 eV. This solution will be further tested by a more 
precise determination of the antineutrino spectrum in KamLAND. 

In any case, our data analysis permit us to considerably limit the size of the a coefficients 
(see Eq. (I26J1 ) and from that to derive a limit on the product of the effective neutrino-scalar 
and matter-scalar Yukawa couplings depending on the mass of the scalar field (Eq. (|27|l). In 
particular, for neutrino-scalar couplings \ u > 3.0 x 10~ 7 (ms/10~ 7 eV) 2 our analysis of solar 
and KamLAND data yields a more restrictive constraint on the matter-scalar couplings than 
gravitational ISL tests. 

These scenarios will be further tested by the precise determination of the energv depen- 
dence of the survival probability of solar neutrinos, in particular for low energies j35j. 
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